arXiv:1509.07970v2 [gr-qc] 30 Dec 2015 


Axially symmetric dissipative fluids in the quasi—static approximation 

L. HerrergQ 

Escuela de Fisica, Facultad de Ciencias, Universidad Central de Venezuela, Caracas, 

Venezuela and Instituto Universitario de Fisica Fundamental y Matemdticas, Universidad de Salamanca, Salamanca, Spain 

A. Di Priscclil 

Fscuela de Fisica, Facultad de Ciencias, Universidad Central de Venezuela, Caracas, Venezuela 

J. Ospinc0 

Departamento de Matemdtica Aplicada and Instituto Universitario de Fisica 
Fundamental y Matemdticas, Universidad de Salamanca, Salamanca, Spain 

J.Carolil 

Departament de Fisica, Universitat Hies Balears, F-07122 Palma de Mallorca, Spain 

Using a framework based on the 1 + 3 formalism we carry out a study on axially and reflection 
symmetric dissipative fluids, in the qnasi-static regime. We first derive a set of invariantly defined 
“velocities”, which allow for an inambiguous definition of the quasi-static approximation. Next we 
rewrite all the relevant equations in this aproximation and extract all the possible, physically rele¬ 
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I. INTRODUCTION 

In the study of self-gravitating fluids we may consider 
three different possible regimes of evolution, namely: the 
static, the quasi-static and the dynamic. 

In the static case, the spacetime admits a timelike, 
hypersurface orthogonal. Killing vector. Thus, a coordi¬ 
nate system can always be choosen, such that all metric 
and physical variables are independent on the time like 
coordinate. The static case, for axially and reflection 
symmetric spacetimes, was studied in [y. 

Next we have the full dynamic case where the system 
is considered to be out of equilibrium (thermal and dy¬ 
namic), the general formalism to analyze this situation, 
for axially and reflection symmetric spacetimes was de- 
velopped in [2| using a framework based on the 1 + 3 
formalism [3l-l6|. 

In between the two regimes described above, we have 
the quasi-static evolution. 

As is well known, in this regime the system is assumed 
to evolve, although sufficiently slow, so that it can be 
considered to be in equilibrium at each moment (Eqs. 
(IMlinil) are satisfied, in the corresponding static case). 
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This means that the system changes slowly, on a time 
scale that is very long compared to the typical time in 
which the fluid reacts to a slight perturbation of hydro¬ 
static equilibrium. This typical time scale is called hy¬ 
drostatic time scale B-i (sometimes this time scale is 
also referred to as dynamical time scale, e.g. i). Thus, 
in this regime the system is always very close to hydro¬ 
static equilibrium and its evolution may be regarded as 
a sequence of equilibrium models. 

Briefly speaking, all the relevant characteristic times 
of the system under consideration, should be much larger 
than the hydrostatic time. 

This assumption is very sensible because the hydro¬ 
static time scale is very small for many phases of the life 
of the star It is of the order of 27 minutes for the 
Sun, 4.5 seconds for a white dwarf and 10“'^ seconds for 
a neutron star of one solar mass and 10 Km radius. It 
is well known that any of the stellar configurations men¬ 
tioned above, generally (but not always), changes on a 
time scale that is very long compared to their respective 
hydrostatic time scales. 

In the spherically symmetric case there exist several 
studies on the behaviour of fluid distributions in the 
quasi-satic regime (see (lol - [l^ and references therein). 

It is our purjTOse here, to make use of the framework 
developped in [^, to carry out a study of axially and 
reflection symmetric fluids in the quasi-static regime. 

For doing that we shall need to introduce different 
invariantly defined “velocities”, in terms of which the 
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quasi-static approximation (QSA) is expressed. 

As we shall see, the shear and the vorticity of the fluid, 
as well as the dissipative fluxes, may affect the (slow) 
evolution of the configuration, as to produce “splittings” 
within the fluid distribution. 

It will be also shown that in the QSA, the contribu¬ 
tions of the gravitational radiation to the components 
of the super-Poynting vector do not necessarily vanish. 
However, as we shall show below, it appears that if at 
any given time, the magnetic part of the Weyl tensor 
vanishes, then it vanishes at any other time afterwards. 
Thus we should not expect gravitational radiation from 
a physically meaningful system, radiating for a finite pe¬ 
riod of time (in a given time interval) in the QSA. 


S, with components 

/ -J 

i^„ = (0,H,0,0); 0,0,A-^- ,0|, 


A 


( 4 ) 


L“ = 


^ ,0, A =,0), (5) 




5« = (0,0,0,C), 

satisfying the following relations: 


( 6 ) 


= -S°‘Sc, = -1, (7) 


II. THE METRIC AND THE SOURCE: BASIC 
DEFINITIONS AND NOTATION 

We shall consider, axially (and reflection) symmetric 
sources. For such a system the most general line element 
may be written in “Weyl spherical coordinates” as: 


ds^ = -A^df+B'^ {dr^ + r^de^)+G^d(l)^+2Gdedt, ( 1 ) 

where A, B, G, G are positive functions of t, r and 9. We 
number the coordinates = t,x^ = r, x^ = 0,x^ = 4>. 

We shall assume that our source is filled with an 
anisotropic and dissipative fluid. We are concerned with 
either bounded or unbounded configurations. In the 
former case we should further assume that the fluid is 
bounded by a timelike surface E, and junction (Darmois) 
conditions should be imposed there. 

The energy momentum tensor may be written in the 
“canonical” form, as 

Tap = (M + P)yaVp + P daP + Aap + QaVp + QpVa- (2) 

The above is the canonical, algebraic decomposition 
of a second order symmetric tensor with respect to unit 
timelike vector, which has the standard physical mean¬ 
ing when Tap is the energy-momentum tensor describing 
some energy distribution, and the four-velocity as¬ 
signed by certain observer. 

With the above definitions it is clear that p, is the en¬ 
ergy density (the eigenvalue of Tap for eigenvector I/“), 
Qa is the heat flux, whereas P is the isotropic pressure, 
and Aap is the anisotropic tensor. We emphasize that 
we are considering an Eckart frame where fluid elements 
are at rest. 

Since we choose the fluid to be comoving in our coor¬ 
dinates, then 



VaK°‘ = V°‘La = V°‘Sa = K°‘La = K°‘Sa = S°‘La = 0. 


( 8 ) 

The unitary vectors U“, form a canonical or¬ 

thonormal tetrad (say ei“^), such that 


„(0) _ y (1) _ (2) _ r (3) _ Q 


with a = 0, I, 2, 3 (latin indices labeling different vec¬ 
tors of the tetrad). The dual vector tetrad is easily 
computed from the condition 


V{a){b) — 9ape-(a.)^^(b)^ 

where Tj{a){b) denotes the Minkowski spacetime metric 
The anisotropic tensor may be expressed through three 
scalar functions defined as (see i): 

Akl = K'^L^Tap, , (9) 

H/ = {2K°‘K‘^ - L°‘L^ - S°^S^)Tap, (10) 


H/j = {2L°‘L^ - S°‘S^ - K°‘K!^)Tap. (11) 

This specific choice of these scalars is justified by the 
fact, that the relevant equations used to carry out this 
study, become more compact and easier to handle, when 
expressed in terms of them. 

Finally, we may write the heat flux vector in terms of 
two scalar functions: 

= qiK^ + qiiLfj, (12) 

or, in coordinate components (see i) 


q^ = 


qilG 


<?/ 


Aqii 


A^A^B^r^ + G2 ’ B ’ ^A'^B^r^ + G^ ’ 


Next, let us introduce the unit, spacelike vectors K,L, 


(13) 
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I y/A^B^r^ + G^qu \ 

0,5^7,---,0 . (14) 


The above scalars may be written in terms of the metric 
functions and their derivatives as (see 0): 


2a I + an = 


3 IB 
A I B 


C 

C 


( 21 ) 


Of course, all the above quantities depend, in general, on 

t,r,e. 

The kinematical variables play an important role in 
the description of a self-gravitating fluid. Here, besides 
the four acceleration, the expansion scalar and the shear 
tensor, we have a component of vorticity. 

Thus we obtain (see 0)- 

For the four acceleration 




(15) 


with 


ai = 


For the expansion scalar 

0 = 

AB^ 


r^A^B'^ + G2 
G2 


'B C 
^B^C 


BAG 
\b~A~^G 


C 

C 


Next, the shear tensor 




where 


/ 7 « = 5 “ + 

may be defined through two scalar functions, as: 


^aj3 — 3^^^^ 


+ g(2cr77 + ai){LaLfj — —hap). 


(16) 


(17) 


(18) 


(19) 


( 20 ) 


2an + ai = 


A'^B'^r'^ + G2 




( 22 ) 


^f_A G_G 

yA~^G~c^ 

where the dot and the prime denote derivatives with re¬ 
spect to t and r respectively. Once again, this specific 
choice of scalars, is justified by the very conspicuous way, 
in which they appear in the relevant equations (see the 
Appendix in 0 )- 

Finally, the vorticity may be described, either by the 
vorticity vector a;“, or the vorticity tensor defined 
as: 


A' _ 

A 

A,61 




AB’ 

Vr^A^B^ + G2 

A 

A2 




to. 


= I i , (23) 


where 0^^ = V^a-,p\ + 0[aVa]) and rjap^n, denote the Levi- 
Civita tensor; we find a single component different from 
zero, producing: 


f^a/3 — f^(TcKAr^ T^TFq,), 


and 


with the scalar function H given by 

ri(G' 2 A'\ 
G A ! 


n = 


2By/A‘^B^r^ +G^' 


(24) 

(25) 

(26) 


Observe that from (ESI) and regularity conditions at 
the centre, it follows that: G = 0 fl = 0. 

Let us now introduce the electric {Ea/s) and magnetic 
(Hafi) parts of the Weyl tensor ( GajS-ys), defined as usual 
by 


Eo.^ = G^^psV''V\ 




P8 


•^PyvyS 


(27) 


The electric part of the Weyl tensor has only three in¬ 
dependent non-vanishing components, whereas only two 
components define the magnetic part. Thus we may also 
write 


EaP = 2 + £ll){EaKp — -hafi) + -I- £l){LaLp — -hap) + SKL{KaLp + KpLa), 


(28) 


















and 


and 
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Ho.p = H,{S^Kp + SpK^) + + SpL^). (29) 

Also, for the Riemann tensor we may define three ten¬ 
sors Yap, Xap and Zap as 

Yap = Ra.psV'V^, (30) 

^a/3 = \ria.^''KpP5V''V\ (31) 


Zap = leaepRsp^^V^, (32) 

where = Irj^pusRap”■ 

From the above tensor, we may define the super- 
Poynting vector by 

Pc = eapj {Y^ZP^ - X]Z^P) , (33) 

which, in our case can be written as: 

Pa = PiKa + Pi I La, (34) 

with 


Pi = —^{2£ii + £i) + 2 Hi£kl + 

Pii = - ~ ‘^^2£kl + 


In the theory of the super-Poynting vector, a state of 
gravitational radiation is associated to a non-vanishing 
component of the latter (see [i^il). This is in 
agreement with the established link between the super- 
Poynting vector and the news functions 2l| , in the con¬ 
text of the Bondi-Sachs approach [23|. Furthermore, 
as it was shown in [^ . there is always a non-vanishing 
component of P^, on the plane orthogonal to a unit vec¬ 
tor along which there is a non-vanishing component of 
vorticity (the 9-r plane). Inversely, P^ vanishes along 
the (((-direction since there are no motions along this lat¬ 
ter direction, because of the reflection symmetry. 

Therefore we can identify three different contributions 
in (1351) . On the one hand we have contributions from 
the heat transport process. These are in principle inde¬ 
pendent of the magnetic part of the Weyl tensor, which 
explains why they remain in the spherically symmetric 
limit. 

On the other hand we have contributions from the 
magnetic part of the Weyl tensor. It is reasonable to 
associate these with gravitational radiation. These are 
of two kinds. On the one hand contributions associated 
with the propagation of gravitational radiation within 
the fluid, and on the other, contributions of the flow of 
super-energy associated with the vorticity on the plane 
orthogonal to the direction of propagation of the radia¬ 
tion. Both contributions are intertwined, and it appears 
to be impossible to disentangle them through two inde¬ 
pendent scalars. 

It is worth noticing that the factors multiplying the 
H terms in (1351) . are £i,£ii,£kl, implying that purely 
magnetic or purely electric sources, do not produce grav¬ 
itational radiation. This is consistent with the result ob¬ 
tained in vacuum for the Bondi metric [2^, stating that 


327r^(7/ 

3 

327r^g/7 


[3(/r + P) -f 11/] -I- 32Tr'^ qiiUkl, 
[3(^ -I- P) -I-11//] -I- 3271^qiIIkl- 


(35) 


purely electric Bondi metrics are static, whereas purely 
magnetic ones, are just Minkowski. 


III. THE KINEMATICS 

The following discussion heavily relies on the kinematic 
quantities characterizing the motion of a medium pre¬ 
sented in [ 2 ^, with slight changes in notation. 

In Gaussian coordinates, the position of each particle 
may be given as 


:r“=x“(y“,s), 


(36) 


where s is the proper time along the world line of the 
particle, and y°‘ (with a running from 1 to 3) is the posi¬ 
tion of the particle on a three-dimensional hypersurface 
(say S). Then for the unit vector tangent to the world 
line (the four-velocity) we have 


and observe that 


ya 


dx°^ 


d ^ _ 1 ^ 


(37) 


(38) 


Next, for an infinitesimal variation of the world line we 
have 

(39) 

from which it follows 


Dt{5x°‘) = V.;^p5x^. 


(40) 
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We can define the position vector of the particle y°‘ + 
5y°‘ relative to the particle 2 /“ on S, as 

5±x°‘ = . (41) 

Then the relative velocity between these two particles, is 

= h^pDriS^x^), (42) 

and considering (|40|) and (ITO it follows that 

= VpjLX^. (43) 

Now, the infinitesimal distance between two neighbor¬ 
ing points on E is 

Sl^ = gayS^x^6^x‘^, (44) 

then 

SlDriSl) = gai3S±x^DTi6±x°'), (45) 

or, by using (HOI) and (H^ . 

SIDt{SI) = Va-isS^xf^S^x^^. (46) 

Then, taking into consideration the expression for the 
irreducible components of a timelike vector 

ka;/3 — 4- ~0/ia/3 T 

and, introducing the spacelike triad 

ef,)^(if“,L“,5“) = (^)(,), (48) 

we obtain 

+ ^hay + fla/?), (49) 

where triad indices i, j run from 1 to 3. 

From ((49l) we can define the following “velocities”: 


1^(1) = (50) 

V(2) = L'^L^ {Cafi + -^&hap + fla/s), (51) 

^(3) = S°^S^{aap -f —Qhap + f^a/?)) (52) 

^(1,2) = K°^ {actP + -^Qhap + (53) 

^(1,3) = K‘^ {(Tap + + f^ct/s); (54) 


which become, using (Uni) and (l24l) 

^(i) = ^(fT/ + 0), F(2) = i(a// + 0), (55) 

V( 3 ) = -^(0 - o-/ - an), V( 1 _ 2 ) = V( 1 ^ 3 ) = 0, (56) 

satisfying 


1^(1) + ^( 2 ) + 1^(3) — 0- (57) 

It is worth noticing that the quantities defined above, 
describe variations of SI, with respect to proper time, 
(projected on different pairs of triad vectors), divided by 
SI. Accordingly, it would be perhaps more appropriate 
to call these quantities “velocity contrast” or “specific 
velocities”. However for simplicity we shall refer to them 
just as “velocities”. 

On the other hand, the geometrical and physical mean¬ 
ing of such quantities, becomes evident from (ISOt - dMl) . 


IV. THE QUASI-STATIC REGIME 

Let us now translate the QSA defined in the first sec¬ 
tion, into conditions to the different definitions of veloc¬ 
ities given above, and kinematical variables. 

The fact that any characteristic time scale of the prob¬ 
lem under consideration must be much smaller than the 
hydrostatic time, implies that: 

• The “velocity” functions V(i)_( 2 )^( 3 ) and V(i, 2 ) de¬ 
fined in (ISOL (SSI are small quantities (say of order 
0 (e), where e << 1 ). 

• We shall neglect all quantities of order O(e^) and 
higher. 

• From (15511551) it follows that tr/ 77 , 0, Q, are of order 
0 (e). 

• From (ITT)) . (1TT1) . (imi . (1T51) it follows that B, C and 
G are of order 0(e). 

• From (dll), (dH), dUl) it follows then, that up to the 
order 0(e) we have that aj = an = a. And up to 
the same order. 


0 -I- (t = 


3B 

AB 


0 - 2 d = 


3C 

'AC' 


(58) 


Next, we have also to assume that the relaxation time 
in the transport equation (Eq.(57) in Q), must be ne¬ 
glected. Indeed, the relaxation time is the time taken by 
the system to return spontaneously to the steady state 
(whether of thermodynamic equilibrium or not) after it 
has been suddenly removed from it. But as it follows 
from the very nature of the QSA, all processes evolve on 
time scales which are much larger than the time scale on 
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which transient phenomena take place, implying that we 
are assuming the heat flux vector to describe a steady 
heat flow. 

Then, neglecting the relaxation time in the trans¬ 
port equation, we obtain the following two equations 
(Eqs.(58,59) in Q, with t = 0) 


qii = - 


A 


Gf + A^Tg 
ABr 


AT a 


II 


(59) 


and 


qi = -^{r + BTai). (60) 

JD 

Therefore, in the quasistatic regime, we obtain from 
the above equations, using the fact that T is of order 
0 (e), and the conditions of thermal equilibrium, 

{TAY = {TA),0 = 0, (61) 

that A is of order 0(e), whereas G is of order O(e^), which 
in turn implies that up to the order 0(e): 


From (1261) 1 it follows at once that in our regime (l is of 
order O(e^) (a result that can also be obtained from B5 
ini). 

From (B6, B7) in i, in the QSA, we have respectively 


Ab' _ a _ a I 

W Br Br \ A ^ G J 


a 

W 


aG' 

BC 


= 87rg/, 

(63) 


2 ^ n' n (2A' G' 


BCr 


Sirqu. 


(64) 

from which it follows that dissipative fluxes are also of 
order 0(e). 

Thus sumarizying all the consequences derived so far 
from the QSA we have: 


• (7, 0,n, B, C, A , qj, qjj, ill, iiii and G are of 
order 0(e). 


• is of order O(e^). 


ai = 


A' ^ 
AB’ 


an 


ABr 


(62) Next, from A6 in i, in the QSA. 
I 


^ + (/^ + F’)0 H—^ (2(T/ -|- ajj) H—^ (2cr// -|- aj) + {qufi + ‘^qiian 


"C" 

1 (BrY] 

qii 

(B,e 


A 

Br 

^ Br 

KB - 

^ G J 


(65) 


it follows that fi is of order 0(e). 


J 


From A7 in i, we obtain the following two equations 
(in the QSA) 


B\ 3 J Br 


G 

Akl,0 + -^Akl 


A. 

'3B 


2G' 

~C~ 


(Br)' 

Br 


Al 

3B 


A 


(Br)' 


Br 


+ I M 

Akl 

Br 


P 

2B,e 

B 


j a/ -I- 11x^0// 


G 


,+^ = 0 
G A 


( 66 ) 


and 

J 


1 

Br 



-b 


A 

A 2 



Ai f B^g CA An /B^g 
^3Br V B ^ G J ^ 3Br \ B ^ 


H + p -I— an -|- Anrai 

‘2‘Gfi \ Akl 

~C~ J ^ “A 


C" „ Br ' 

-h 2 - - - 

G Br 


^ = 0 
A 


( 67 ) 
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Since the hydrostatic equilibrium condition holds at any 
time, the corresponding hydrostatic equilibrium equa¬ 
tions (Eqs. (21,22) in [2]) must be satisfied. Then, from 
the two equations above, we obtain respectively 

nifL~0(e); qi-0{e^); B ^ 0{e^) C^Oie^). 

( 68 ) 

and 

ri77~0(e); g77~0(e2) P« 0(e), (69) 

where the fact has been used that in the QSA, P, 11/, 11// 
contain (besides terms including spatial derivatives of the 
metric tensor), terms with B and C. 


From (( 68 ll and (l58l) it follows at once that 

aKiQ{e^)- 0ft!O(e2). (70) 

Let us now turn to (l63ll . which, using (l55ll and dSi, 
may be written as: 

2V' = a [In (dC)]' -f ^ [in (nA^C)] ^ + SnqjB, (71) 

with V = V(i) = V( 2 )- 

After integration we obtain 

V = Vs-- / {d [ln(dC')]'-|-— [in(nA^C')] g+Snq/Bjdr, 

(72) 


^(3) ~ Idsly; ~ — 


(3)E 



r ,cA 

/ 

In(^) 

Jr 

a 


+ — [in (flA^C)] g -I- SnqjB [ dr, 


(73) 


where the boundary surface of the source is defined by 
the equation r = rs and the fact that V( 3 ) = V — d, has 
been used. 


In a similar way we may writte down (16411 as 


2 V 0 = d [In (dC)] 0 — i^r 


In 



+ 8 TrqiiBr, (74) 


producing 


J 


^(3) 




where now the boundary surface equation is given by 

9 = 9s. 

Let us focus on the expressions above dza), dza, dza, 
(1761) . If we assume the fluid to be irrotational, shear- 
free and dissipationless, then the sign of V and 1 /( 3 ) is 
the same as the sign of Vs and V( 3 )i;, for any fluid el¬ 
ement within the distribution. However, the presence 
of any of the factors above (vorticity, shear, heat flux), 
may lead to a situation, where either velocity changes 
of sign within the fluid distribution, with respect to its 
sign on the boundary surface. In other words, it may 
happen that some inner regions move in one direction 
whereas the outer ones move in the opposite direction. 
Such “splittings” of the configuration, have already been 
reported for the spherically symmetic case (see [l2| , [III , 


V = Vs-^ ( "" \d[\n{dC)] g-rn [in '-f 87 rg//Pr} 


(75) 


or 


rH [in (A^CH)]^ -|- SirqjjBr > d9, 


(76) 


I- 

[i3[i dll, and references therein). Here the picture 
is more involved than in the spherically symmetric case 
due to the possibility of the splitting to occur, along two 
orthogonal directions. 

Next, Eqs. (B 8 , B9) in Q, read in the QSA as 




I 


2Br 


<1,9 


dC^i 

C 


(77) 


r. 1 

H2 = —flail — „ „ fl,e — 


2Br 


C 


2B 


a V 


~C 


[78) 

implying that the magnetic part of the Weyl tensor is of 
order 0 (e). 


d9, 


















It is worth noticing that in the vorticity-free case( G = 
H = 0), it follows from (1771) and (1751) that 






2BrC ’ 


H2 


2BC ■ 


(79) 


Then, from (1551) and (IMl) . using (1551) . we find (always 
assuming = 0) 


2 



{^cy 

c 


+ SttBqi, 


(80) 


and 


2 



C 


+ S-KBrqij. 


(81) 


The combination of the two equations above with (1751) 
produces 


Hi 


1 

Br 



+ 47rq//, 
.9 


(82) 


V. CONCLUSIONS 

Wew have provided a general framework for describing 
the evolution of axially symmetric dissipative fluids in the 
QSA. 

The role played by the vorticity, the shear and the dis¬ 
sipative flux is clearly brought out, through the expres¬ 
sions dZlI), (1751) and dZll), dZil). Such expressions show 
how the fluid distribution may split, under the effects of 
the factors mentioned above, leading to a variety of very 
different structures. 

Finally it is worth mentioning that in the QSA the 
magnetic part of the Weyl tensor does not necessarily 
vanish (though it is of order 0(e)), thereby implying 
that the “gravitational” part of the super-Poynting vec¬ 
tor does not vanish either, meaning that gravitational 
radiation is not incompatible with the QSA. 

However, from (1771) and (1751) (alternatively see (IA.8I) 
and (IA.9I) 1. it follows at once that Hi Ki H 2 ~ O(e^) (at 
least), and therefore are neglectable in the QSA. This in 
turn implies that if the magnetic part of the Weyl tensor 
vanishes at any given time, it will do so for any time 
afterwards. In other words, no state of radiation for a 
finite period of time is expected in the QSA. This result 
is in agreement with the one obtained by Bondi [s^, for 
a more restricted case. However besides the “inductive” 
transfer of energy, mentioned by Bondi, we also have 
here, the transfer carried on by the dissipative flux. 


and 

Thus from (1791) it follows at once that, in the QSA, 
for the vorticity-free case, the vanishing of the shear is 
a necessary and sufficient condition for the fluid to be 
purely electric. This result is somehow complementary to 
the one obtained in the general dynamic case [30l| , which 
states that, for a dissipative and anisotropic, shear-free 
fluid, the vanishing vorticity, is a necessary and sufficient 
condition for the magnetic part of the Weyl tensor to 
vanish, this last result in turn provides a generalization 
of a similar result for perfect fluids, obtained in [sTI - fs^ . 
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Appendix: The remaining equations in the QSA 

Below, for the benefit of the reader, we shall write the 
equations not used explicitly in the text, specialized for 
the QSA, from the framework developped in 0. Thus, 
Eqs. (B10-B18) in @], read respectively 
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